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To clarify the instability of the ferrimagnetism which is the fundamental magnetism of ferrite, numerical- 
diagonalization study is carried out for the two-dimensional S = 1/2 Heisenberg antiferromagnet with frus¬ 
tration. We find that the ferrimagnetic ground state has the spontaneous magnetization in small frustration; 
due to a frustrating interaction above a specific strength, the spontaneous magnetization discontinuously 
vanishes so that the ferrimagnetic state appears only under some magnetic fields. We also find that, when 
the interaction is increased further, the ferrimagnetism disappears even under magnetic field. 


Ferrite is a magnetic material that is indispensable in modern society. It is be¬ 
cause this material is used in various industrial products including motors, generators, 
speakers, powder for magnetic recording, and magnetic heads etc. It is widely known 
that fundamental magnetism of the ferrite is ferrimagnetism. The ferrimagnetism is 
an important phenomenon that has both ferromagnetic nature and antiferromagnetic 
nature at the same time. The occurrence of ferrimagnetism is understood as a mathe¬ 
matical issue within the Marshall-Lieb-Mattis (MLM) theorem^’concerning quantum 
spin systems. A typical case showing ferrimagnetism is when a system includes spins 
of two types that antiferromagnetically interact between two spins of different types in 
each neighboring pair, for example, an [S, s) = (l, 1/2) antiferromagnetic mixed spin 
chain, in which two different spins are arranged alternately in a line and coupled by the 
nearest-neighbor antiferromagnetic interaction. The ferrimagnetic state like the above 
case, in which the spontaneous magnetization is fixed to be a simple fraction of the 
saturated magnetization determined by the number of up spins and that of down spins 
in the state, is called the Lieb-Mattis (LM) type ferrimagnetism. Another example of 
ferrimagnetism is a system including single-type spins that are more than one in a unit 
cell, although the ferrimagnetism can appear even in a frustrating system including 
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Fig. 1. (Color) Network of antiferromagnetic interactions studied in this paper. The black and red 
bonds represent Ji and J 2 interactions. Green squares denote finite-size clusters of 24 and 30 sites in 
(a) and (b), respectively. Note that the two-dimensional network composed only of the black bonds is 
called the Lieb lattice. 


only a single spin within a nnit cell.^’®^ 

The antiferromagnet on the Lieb lattice illustrated in Fig. 1 corresponds the second 
case, in which there are three spins in a unit cell. The MLM theorem holds in the Lieb- 
lattice antiferro magnet. If antiferromagnetic interactions are added to this Lieb lattice 
so that magnetic frustrations occur, however, the MLM theorem no longer holds. In this 
situation, the ferrimagnetic state is expected to become unstable. The problem of how 
the ferrimagnetism collapses owing to such frustrating antiferromagnetic interactions is 
an important issue to understand the ferrimagnetism well and to make ferrimagnetic 
materials more useful in various products. This problem was studied in the S' = 1/2 
Heisenberg antiferromagnet on the spatially anisotropic kagome lattice,®’where the 
existence of an intermediate phase with weak spontaneous magnetization is clarihed 
between the LM type ferrimagnetic phase and the nonmagnetic phase including the 
isotropic kagome-lattice antiferromagnet. We are then faced with a question: is there 
any other different behavior of the collapse of the ferrimagnetism? 

Under circumstances, the purpose of this study is to demonstrate the existence of 
a different behavior of collapsing ferrimagnetism in the case of an S' = 1/2 Heisenberg 
antiferromagnet on the lattice shown in Fig. 1 to answer the above question. When the 
antiferromagnetic interactions denoted by the red bonds vanish, the system is unfrus¬ 
trated and thus it certainly shows ferrimagnetism in the ground state. In this study, we 
examine the case when the red-bond interactions are switched on. 
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The model Hamiltonian examined in this study is given by "H = 'Ho + 'Hzeeman, where 
Ho = 


Hzeeman 

Here Si denotes an S = 1/2 spin operator at site i. Sublattices a, a', and (3 and the 
network of antiferromagnetic interactions Ji and J 2 are depicted in Fig. 1. Here, we 
consider the case of isotropic interactions. The system size is denoted by N^. Energies 
are measured in units of Ji; thus, we take Ji = 1 hereafter. We examine the properties 
of this model in the range of J 2 /J 1 > 0. Note that, in the case of J 2 = 0, sublattices a 
and a' are combined into a single sublattice; the system satishes the above conditions 
of the MLM theorem. Thus, ferrimagnetism of the LM type is exactly realized in this 
case. In the limit of J 2 / J\ —)■ cxd, on the other hand, the lattice of the system is reduced 
to a trivial system composed of isolated S' = 1/2 spins and isolated dimers of two spins. 
Its ground state is clearly different from the state of the LM-type ferrimagnetism in the 
case of J 2 = 0. One thus hnds that while J 2 becomes larger, the ground state of this 
system will change from the ferrimagnetic one in the case of J 2 = 0 to another state, 
which we survey here. 

Next, we discuss the method we use here, which is numerical diagonalization based 
on the Lanczos algorithm.It is known that this method is nonbiased beyond any 
approximations and reliable for many-body problems, which are not only localized spin 
systems such as the Heisenberg modeh^’ treated in th present study but also strongly 
correlated electron systems including the Hubbard modeh^^^®^ and the t-J model. 

A disadvantage of this method is that the available system sizes are limited to being 
small. Actually, the available sizes in this method are much smaller than those of the 
quantum Monte Carlo simulation^®’and the density matrix renormalization group 
calculation;^^^ however, it is difficult to apply both methods to a two-dimensional (2D) 
frustrated system like the present model. This disadvantage comes from the fact that 
the dimension of the matrix grows exponentially with respect to the system size. In 
this study, we treat the hnite-size clusters depicted in Fig. 1 when the system sizes 
are W = 24 and 30 under the periodic boundary condition. Note that each of these 
clusters forms a regular square although cluster (b) is tilted from any directions along 
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Fig. 2. (Color) Results for J 2 / Ji = 0.55. Lowest energy in each subspace of M for the system of 
Ns = 30 is shown in panel (a). The magnetization process is depicted in panel (b); red and black lines 
represent results for Ns = 24 and 30, respectively. 


interaction bonds. 

We calculate the lowest energy of T-Lq in the subspace characterized by Yhj = Mhj 
numerical diagonalizations based on the Lanczos algorithm and/or the Householder al¬ 
gorithm. The energy is represented by E{Ns, M), where M takes every integer up to the 
saturation value Mg (= SNg). We here use the normalized magnetization m = M/Mg. 
Some of Lanczos diagonalizations have been carried out using the MPI-parallelized 
code, which was originally developed in the study of Haldane gaps.^^^ Note here that 
our program was effectively used in large-scale parallelized calculations. 

To obtain the magnetization process for a hnite-size system, one hnds the magneti¬ 
zation increase from M to M + 1 aX the held 

h = E{Ng,M + l)-E{Ng,M), (3) 

under the condition that the lowest-energy state with the magnetization M and that 
with M + 1 become the ground state in specihc magnetic helds. Note here that it often 
happens that the lowest-energy state with the magnetization M does not become the 
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Fig. 3. (Color) Results for J 2 /J 1 = 0.64. Lowest energy in each subspace of M for the system of 
Ns = 30 is shown in panel (a). The magnetization process is depicted in panel (b); red and black lines 
represent results for TVs = 24 and 30, respectively. Main panel is a zoomed-in view of its inset with a 
wide range. The broken lines represent the results before the Maxwell construction is carried out. 


ground state in any field. The magnetization process in this case is determined around 
the magnetization M by the Maxwell construction.^®’ 

Now, we observe the case of J 2 IJ\ = 0.55; results are shown in Fig. 2. Figure 2(a) 
depicts the lowest energy level in the subspace belonging to M for iVs = 30. The levels 
for M = 0 to M = 5 are identical within the numerical accuracy. For M > 5, the 
energies increase with M. This behavior indicates that the spontaneous magnetization 
is M = 5. In Fig. 2(b), we draw the magnetization process determined by eq. (3) in the 
full range from the negative to the positive saturations. The spontaneous magnetization 
m = 1/3 appears and the state at m = 1/3 shows the plateau with a large width. It is 
observed that, above m = 1/3, the magnetization grows continuously. These behaviors 
are common with those of the LM ferrimagnetism at the unfrustrated case of J 2 = 0. 

Next, let us examine the case of J2/J1 = 0.64; results are shown in Fig. 3. The 
M dependence of the lowest energy belonging to M is different in M < 3 from the 
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case of J 2 /J 1 = 0.55. This difference affects with the disappearance of the spontaneous 
magnetization, which is shown in Fig. 3(b). This discontinuous disappearance occurs 
at J 2 IJ\ ~ 0.59 for iVg = 24 and at J 2 IJ 1 ~ 0.63 for = 30. An important point 
is that an intermediate state with smaller but nonzero spontaneous magnetizations is 
absent between the m = 1/3 state and the nonmagnetic state. This behavior is clearly 
different from the presence of such an intermediate state in the spatially anisotropic 
kagome lattice.®’We speculate that this difference comes from the point that the 
competing interaction in the present model has a strong quantum nature localized 
at pairs of dimerized spins. The discovery of the future third case of the collapsing 
ferrimagnetism would contribute to conhrm our speculation. Note also that the plateau 
at m = 1/3 shows a large width. This suggests that the ferrimagnetic state is realized 
if external magnetic helds are added. 

To examine the properties of the m = 1/3 states in a more detailed way, we evaluate 
the local magnetization dehned as 

“1m = ^ ('i) 

where ^ takes a, a' and jS. Here, the symbol {O) denotes the expectation value of the 
operator O with respect to the lowest-energy state within the subspace with a hxed 
M of interest. Recall here that the case of interest in this paper is M = M^/S. Here 

denotes the number of ^ sites. Results are shown in Fig. 4. In the region of small 
J2/J1, a and a' spins are up and (3 spin is down, although each of magnetizations is 
slightly deviated from the full moment due to a quantum effect. This spin arrangement 
is a typical behavior of ferrimagnetism. On the other hand, in the region of large J 2 / Ji, 
the magnetizations at a and a' spins are vanishing and (3 spin shows almost a full 
moment up. This marked change in the local magnetizations occurs at J 2 /J 1 ~ 1.38 for 
As = 24 and at J2/J1 ~ 1.40 for Ag = 30, which suggests the occurrence of the phase 
transition around at J 2 IJ\ ~ 1.4. Therefore one hnds that, for J 2 /J 1 larger than this 
transition point, the ferrimagnetic state cannot be realized even under magnetic helds. 
It is unfortunately difficult to determine the transition point in the thermodynamic limit 
precisely only from the present two samples of small clusters. For the determination, 
calculations of larger clusters are required in future studies. Note here that similar 
observations of the local magnetizations were reported in Refs. 28 and 29, which treated 
the Heisenberg antiferromagnet on the Cairo-pentagon lattice,a 2D network obtained 
by the tiling of single-kind inequilateral pentagons. The same behavior of ttilm is also 
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observed when the kagome-lattice antiferromagnet^^^^'^^ is distorted in the \/3 x y/S 
type.^^’^®^ The relationship between these models should be examined in future studies. 
Note also that, in the present model, the change around the transition point seems 
continuous irrespective of whether the system size is Ng = 24 or 30. This aspect is 
different from the observation in the Cairo-pentagon-lattice antiferromagnet,^®’^®^ where 
the change around the transition point seems continuous for Ng = 24 but discontinuous 
for Ng = 30. We speculate that whether the change is continuous or discontinuous in 
hnite-size data is related to whether the number of unit cells in hnite-size clusters is an 
even integer or an odd integer. To confirm this speculation, further investigations are 
required in future. It will be an unresolved question whether the transition is continuous 
or discontinuous in the thermodynamic limit. Figure 5 depicts the magnetization process 
at J2/J1 ~ 1.39. No jumps seem to appear in the process at J2/J1 corresponding to 
the transition point. It is unclear whether the width at m = 1/3 survives or vanishes 
although this m = 1/3 width at J2/J1 ~ 1-39 is smaller than those in Figs. 2(b) and 
3(b). Future studies would clarify how the magnetization process behaves in the vicinity 
of the transition point. 

In summary, we have investigated how the ferrimagnetic state of the S = 1/2 
Heisenberg antiferromagnet on the 2D lattice collapses owing to magnetic frustration by 
numerical-diagonalization method. We capture a discontinuous vanishing of the sponta¬ 
neous magnetization without intermediate phase showing spontaneous magnetizations 
that are smaller than that of the Lieb-Mattis ferrimagnetic state when a frustrating 
interaction is increased. We also observe the disappearance of the ferrimagnetic state 
under magnetic helds for even larger interaction showing frustration. It is known that 
organic molecular magnets can realize ferrimagnetism.^®’'^^^ Since variety of lattice struc¬ 
ture leading to an interaction network is available in such organic molecular magnets, the 
experimental conhrmation might be done in these magnets more easily than metallic- 
element compounds. Further studies concerning instability of the ferrimagnetism would 

contribute much for our development of more stable ferrimagnetic materials. 
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Fig. 4. (Color) Behavior of local magnetizations vs. the ratio of interactions J\ together with a 
zoomed-in view near the transition point in inset. Closed circles and closed diamonds denote results 
for a and f3 for Ns = 24, respectively. Results for a' for Ns = 24 are identical those for a within the 
numerical accuracy because a and a' are symmetric in the Ns = 24 cluster. Open circle, open 
triangle, and open squares represent results for a, a', and j3 for Ng = 30, respectively. Due to the 
tilting for Ng = 30, a and a' are not symmetric, although results of a and a' for Ns = 30 are slightly 
different but very similar. To avoid invisibility from overlapping of symbols, results of a' are shown 
only in inset. 
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Fig. 5. (Color) Magnetization process for J 2 /J 1 = 1.39. Red and black lines represent results for 
Ns = 24 and 30, respectively. 
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